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ABSTRACT 

ASCA observations of NGC 4636 and a southern region have revealed 
extended X-ray emission to a radius of about 300 kpc from the galaxy. 
The symmetric nature of the observed surface brightness around NGC 4636 
indicates its association to this galaxy rather than to the Virgo cluster. Model 
independent estimation of the gravitational mass profile shows a flattening 
at a radius of 20 ~ 35 kpc, where the total mass reaches ~ 6 x 1O 11 M and 
a mass-to-light ratio of 23. The mass still increases to larger radii, reaching 
9 x 10 12 M Q with a mass-to-light ratio of 300 at ~ 300 kpc from NGC 4636. 
These features suggest presence of a galaxy group surrounding NGC 4636. If 
such optically dark groups are common among X-ray bright ellipticals, it would 
explain the very large scatter in their X-ray luminosities with similar optical 
luminosities. 

Subject headings: galaxies:ISM — X-rays:galaxies — galaxydndividual 
(NGC 4636) 
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1. Introduction 

Spatial distribution of gravitational mass around elliptical galaxies has been estimated 



from radial distribution of X-ray emitting hot interstellar medium (ISM) (e.g. Fbrman et 



al. 19851 [lrinchieri et al. 1986|). However, the outermost radius of the ISM distribution has 



so far remained undetermined due to a lack of sensitivity. As a result, the total gravitating 
masses of elliptical galaxies depend heavily on the assumed overall extent of the ISM. High 
sensitivity X-ray observations of bright galaxies are therefore important in determing the 
mass distribution in these systems. 

NGC 4636 is one of the nearest giant elliptical galaxies. We employ a distance of 
17 Mpc ( [TuHy 1988|) . Diffuse thermal X-rays from its ISM make it one of the X-ray 
brightest elliptical galaxies (e.g. [Forman et al. 1985| ; [Trinchieri et al. 1986|) except those 



located at the center of clusters. Although NGC 4636 is located in the Virgo Southern 
Extension (|Nolthenius 1993]) , it is apparently clear of the large-scale X-ray emission from 



hot intra-cluster gas associated with the Virgo cluster ( Takano et al. 1989 ; Bohringer et al 



1994| ). Trinchieri et al. (1994), using ROSAT, detected a largely extended X-ray emission 
up to 18' from NGC 4636, and suggested it to be an extension of the Virgo cluster emission. 
Here, we present new results on the extended X-ray emission around NGC 4636 based on a 
very long observation towards the galaxy center and on an offset observation from ASCA. 
Spectral properties of NGC 4636 have been already reported by Matsushita et al. (1997). 



2. Observation 

We observed NGC 4636 in the 0.5-10 keV X-rays twice with the ASCA observatory 
( I'lanaka et al. 1994| ). The first observation, performed for 40 ksec in 1993 July 22, was 
already analyzed by several authors ( |Awaki et al. 199"4] ; |Matsushita et al. 1994| ; Mushotzky 
|et al. 1994"t ; |Matsumoto et al. 1997| ; [Matsushita et al. 1997|) . The second observation was 
conducted from 1995 December 28 through 1996 January 4, resulting in an exceptionally 
long exposure for 200 ksec. We also employ ASCA data taken on 1994 June 15, in the field 
of NGC 4643, which is just to the south of NGC 4636. 



-4 - 



3. Data Analysis and Results 

3.1. X-ray image and radial brightness profile 

Figure 1 shows background-subtracted 0.7-2.0 keV X-ray images of the NGC 4636 
field and the south-offset field, taken with the GIS (Gas Imaging Spectrometer) ( |Ohashi et 



aL 1996| ; [Makishima et al. 1996|) . This energy band is dominated by the ISM emission. In 



addition to the bright emission peaked at NGC 4636, a fainter emission spreads over the 
entire GIS field of ~ 25' radius, and even to the south-offset field ~ 1° (~ 300 kpc) away. 

The extended emission looks azimuthally symmetric around NGC 4636, without 
significant brightening toward the Virgo center (to the north); along a circle of a radius 
10' centered on NGC 4636, the brightness is constant to within 30%. Furthermore, Ginga 
showed an upper limit on the diffuse X-ray intensity at ~ 2° north of NGC 4636 ( [Takano| 
et al. 19S^| ) to be 10 times lower than that observed in Fig. 1 at the edge of the NGC 4636 



field. We therefore conclude that the extended emission is associated with NGC 4636 rather 
than the Virgo cluster. 

We derived radial profile of X-ray surface brightness, as shown in Fig. 2, by azimuthally 
averaging the intensity in the 0.7-2.0 keV energy range around NGC 4636. The two 
bright background sources in the north of NGC 4636 were excluded, while effects of fainter 
background sources are negligible (< 5%) flTrinchieri et al. 1994j ). The discontinuity of the 
profile at ~ 25', corresponding to the field boundary, is caused by the vignetting and stray 
light of the ASCA X-ray Telescope (XRT; fSerlemitsos et al. 1995| ; [1'akahashi et al. 1995| ). 
As shown in Table 1 (the last column), the fraction of stray light from the NGC 4636 
field to the south offset field is about 47% and 42% at 25' < r < 50' and at 50' < r < 80' 
respectively, where r is the radius measured from NGC 4636. If we limit the region to 
60' < r < 80', the scattered light from the r < 50' and r < 60' field amounts to 64% and 
80% of the detected photons, respectively. As shown later, we conclude from the brightness 
profile in Fig. 2 that the X-ray emission is significant until 60' (300 kpc) from NGC 4636. 

We fit the radial X-ray profile with f3 models convolved with the instrumental 
response of ASCA. The (3 model describes projected X-ray surface brightness S{r) as 
S(r) = So jl + ( r / r coro) 2 } 13 2 , where r is the projected radius, while r core (the core 
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radius), So, and (3 are parameters (e.g. Sarazin 1988). Based on the previous Einstein 
results ( [Forman et al. 1985| ; |Trinchieri et al. 1986Q , we fixed r core to 10" which is practically 



unresolved with ASCA. As shown in Fig. 2, the data within r ~ T (~ 35 kpc) can be 
approximated well by the (3 = 0.5 model, which behaves essentially as oc r~ 2,0 . In contrast, 
the data start deviating from the model beyond R ~ 7'. These properties agree very well 
with the ROSAT radial X-ray profile ( [Trinchieri et al. 1994|) . Clearly, the observed radial 
X-ray profile cannot be fitted with a single /3-model. 

We accordingly fitted the radial profile with a sum of two (3 models, a compact and 
an extended ones, both centered on NGC 4636. Except for the core radius of the compact 
one, the two components were allowed to have free parameters. The radial extent of the 
extended component was examined by changing a cut-off radius, and models with its value 
less than 50' predicts too little flux at r > 50' with the reduced x 2 > 1-75 for 64 dof.. When 
we assumed the cut-off radius to be 60' and 150', the reduced \ 2 decrease to 1.34 and 1.12, 
respectively. This shows that the X-ray emission is significant to 60' from NGC 4636. We 
assumed the cut-off radius to be 150' in the later analysis. As shown in Fig. 2, the double f3 
model with this cut-off radius can fit the observed profile. The (3 values are 0.60 (0.57-0.63) 
and 0.40 (0.38-0.42) for the compact and extended components respectively, and r core of 
the extended component is 6' (4' — 8') or 30 (20-40) kpc. The two (3 components cross over 
at r ~ 5' (~ 25 kpc). Since we can exclude a flat profile for the wider component, it is 
certain that the extended emission is localized in the region surrounding NGC 4636. When 
integrated up to r = 60', the 0.5-4 keV luminosity carried by the extended (3 component 
becomes 8.1 x 10 41 ergs s -1 , which exceeds that of the compact (3 component by a factor of 
5. 



3.2. Temperature profile of the ISM 

In order to study temperature distribution of the gas, we accumulated the GIS data 
in concentric annular regions as summarized in Table 1. Then we can obtain a projected 
(not a three-dimensional) temperature profile. Spectral fit is carried out with the R-S 
model with photoelectric absorption. To represent the contribution from low-mass X-ray 
binaries in NGC 4636 (Matsushita et al. 1994), we added a thermal bremsstrahlung with 
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a temperature fixed at 10 keV. Since GIS data cannot constrain hydrogen column density 
(iVn) accurately, we fixed the value to those determined by the SIS in Matsushita et al. 
(1997). In the same way as in Matsushita et al. (1997), the abundances of heavy elements 
are represented by two parameters, one for the a-elements (O, Ne, Mg, Si, and S) and the 
other for the Fe (Fe and Ni) group, which are denoted as A a and Ap e respectively. At each 
radius we fixed A a to the best fit value derived by Matsushita et al. (1997), and left Ap e 
free. Beyond r > 12', we fixed A a to 0.37 solar as derived for r = 8' — 12' by Matsushita 
et al. (1997), because the field of view of the SIS does not cover this region and no good 
information is available. 

Acceptable fits have been obtained, and the results are summarized in Table 1. The 
spectral fits indicate that the gas is approximately isothermal, but there is a small positive 
temperature gradient, from 0.7 to 0.8 keV, from the center to ~ 10', and beyond which the 
temperature drops to 0.5 keV at the south-offset field (r > 25') (Table 1). The ROSAT 
data show steeper gradient within several arcminutes from the center, probably because of 
the higher angular resolution. 

This "ring-sorted analysis" is however affected by the spectral mixing effect among 
different sky regions, due to the extended point spread function (PSF) of the XRT as well as 
the projection effect. The PSF tends to smear out a radial gradient in the ISM properties. 
Fraction of leakage photons from other annular regions into each ring estimated from the 
observed brightness distribution is about 50% (2nd column from the right in Table 1) for 
r = 5' ~ 50'. 

3.3. Mass estimation from the double-/? model 

The total mass M(R) within a 3-dimensional radius R, assuming a 
hydrostatic equilibrium and a spherical symmetry, is given by (e.g [Sarazin 1988 ), 
M(R) = — ( fjfil + dinfl ) ' wnere m p i s ^e proton mass, k is the Boltzmann constant, 
G is the constant of gravity, and /i ~ 0.63 is the mean particle mass in unit of m p . 

To estimate the radial profile of the total mass, we analytically deproject the two 
/3-model components in Fig. 2 back to 3 dimensions individually, and obtain emissivities of 
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the compact and extended components, e c and e e , respectively. We then calculate the ISM 



density profile as the sum of the two components; n(R) = y{e c (i?) + e e (R)}/A(T, A) , 
with A(T, A) denoting the emissivity function for a temperature T and metal abundance A. 
We assume the two components to have the same temperature. 

Based on the projected temperature and abundance profiles as shown in Table 1, we 
approximate the 3-dimensional profiles for temperature by T(R) = 0.8 — 0.3exp(— R/R^) 
keV with R^ = 10 kpc QTrinchieri et al. 1994j ), and for the abundance profile by 



A(R) = 0.2+0.8 exp(-R/R A ) solar with # A = 40 kpc flMatsushita et al. 1997|) . The observed 
temperature drop beyond 50 kpc is modeled as T(R) = 1.1 — 0.3exp{(_R — R^J/R^} keV 
with i?T a = 100 kpc and i?T b = 400 kpc. These approximation are confirmed to reproduce 
the observed projected profiles very well. 

The derived total mass of NGC 4636 is presented in Fig. 3, together with the ISM 
mass profile (obtained by integrating n) and the stellar mass profile assuming the r 1//4 -law 
(de Vaucoulers 1948) with assuming M/L (mass-to-light ratio) of stars to be 8. We also 
show the result for a constant temperature (0.80 keV) and abundance (1.0 solar) in Fig. 3. 
The total mass M(R) increases roughly as oc R up to ~ 10 kpc, and then it once flattens at 
(3 — 4) x 10 11 Mq, producing a shoulder-like structure at R ~ 20 kpc. Then M(R) starts 
increasing again beyond R ~ 30 kpc. This is essentially the same feature as seen in the 
Fornax cluster fllkebe 1995| ; [Ikebe et al. 1996|) and A1795 (|Xu et al. 1997] ). At R ~ 300 kpc, 



M(R) reaches 9 x 10 12 M Q . 



3.4. Model-independent mass estimation 

The step-like feature in the gravitational mass profile derived in the previous section 
may depend on the particular choice of the double-/? model. We need to examine if other 
mass profiles are consistent with the observed brightness profile. To answer this question, 
we estimate the total mass without assuming the double-/? model. For simplicity, we assume 
an isothermal gas distribution since variation of emissivity with the gas temperature in the 
range 0.5~0.8 keV is less than 20%, and effects of temperature and abundance gradients on 
the total mass were found to be relatively small. From Fig. 2, we can safely assume that 
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beyond 20' and within r < 3', the radial brightness is dominated by a single /? component 
with (3 = 0.4 and /? = 0.6 with r core = 10", respectively. These features are consistent with 
the Einstein and ROSAT results. Therefore, in the regions R < 3' and R > 20', we assume 
that the gravitational mass profile is same as that from the double-/? analysis. Then, our 
task is to find the mass profile which can smoothly connect the two regions. 

We searched for the mass profile in the R — 3' ~ 20' range which can produce the 
observed brightness profile by means of Monte-Carlo simulations. We divided the radius 
into n (n — 3,4,5,6,8,9,10 and 16) equal intervals in a logarithmic scale and assumed 
that d log M/d log R is constant within each step. We randomly picked up a value of 
d log M/d log R between R = 3' and R = 20' for 2 x 10 6 times for each of the n radial 
intervals. In the actual analysis, we extended the mass profile to R < 3' with the same 
d log M/d log it! value at R — 3'. For each trial mass curve, we compare the implied 
brightness profile with the observed one through x 2 fits. We take only acceptable results 
at 90% confidence limits. After a number of trials, 153 acceptable mass curves have been 
obtained. The envelope of these curves, is plotted in Fig. 3, which would represent the 
acceptable range of the gravitational mass from this model independent analysis. We note 
that M(R) for the best fit double-/? model lies within the envelope of model-independent 
estimation. Although the uncertainty of M(R) became larger than that of double-/? 
modeling, variation of M(R) holds some similarity. The growth rate of M(R) drops at 
R ~ 10 — 30 kpc, and starts increasing again in the outer region. 

4. Discussion 

ASCA observations have shown a very extended X-ray emission with a radius of 
at least 60' (300 kpc) surrounding NGC 4636, which is much larger than that detected 
by ROSAT. We studied radial temperature and density distributions of the gas, and 
determined gravitational mass profile with and without assuming a double-/? model. 

According to Fig. 3, the gravitational halo of NGC 4636 appears to terminate once at 
R ~ 10 — 30 kpc, where the total mass reaches several xlO 11 M Q . This inner component 
probably corresponds to the mass associated with the galaxy NGC 4636. Then, the total 
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mass starts increasing again, forming a halo-in-halo structure whose total size is comparable 
to that of a galaxy group. The mass of the whole system and its mass-to-light ratio at 300 
kpc reach ~ 1 x 10 13 M and ~ 300, respectively, which are again comparable to those of 
galaxy groups ( |Mulchaey et al. 1996| ). 



At R > 200 kpc, the X-ray emitting plasma becomes the dominant form of baryons, 
with its mass reaching 5-8% of the total gravitating mass. This value is considerably higher 
than those of individual elliptical galaxies, and again close to those found in galaxy groups 
and poor clusters flMulchaey et al. 1996|) . These features all support the presence of a 
gravitational potential with a size of a galaxy group around NGC 4636. A further support 
is given by the observed abundance decrease in the X-ray emitting plasma flTrinchieri et 
a [. 1994| ; |Mushotzky et al. 1994| ; |Matsushita et al. 19971 ; IMatsushita T997D , from ~ 1 solar 



within ~ 30 kpc, to ~ 0.2 solar beyond ~ 50 kpc to a is typical level of groups and clusters 
of galaxies. 

Nolthenius (1993) identified NGC 4636 as a member of Virgo Southern Extension F 
Cloud. However the 3-dimensional location of NGC 4636 is far offset from the center of the 
member galaxies in the cloud, and their velocity dispersion, 463 km s _1 , is much larger than 
that inferred from the observed gas temperature. The diameter of the cloud, 1.26 Mpc, is 
also much larger than the scale of the X-ray emission. Therefore, the relation between the 
X-ray halo and the galaxy cloud remains unclear, but it is very likely that there is some 
galaxy concentration around NGC 4636. 



Extensive studies of ASCA data have revealed (|Matsushita 19970 that X-ray luminous 
elliptical galaxies preferentially possess large-scale X-ray halos of a few 100 kpc scale, just 
like NGC 4636. Such a galaxy may be regarded as a dominant member of a galaxy group 
even though its evidence is scarce in the optical data. Thus, the X-ray emitting plasma 
may provide a better tracer of the total mass distribution than the light emitting matter. 

We speculate that X-ray luminous elliptical galaxies may commonly possess such 
an extended emission with low surface brightness. X-ray luminosity of elliptical galaxies 
scatter by 2 orders of magnitude for the same optical luminosity, and its cause has long 
been a puzzle. The presence and absence of the extended emission can easily account 
for the large scatter in the total X-ray luminosity, and the low surface brightness of the 
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extended emission, such as in NGC 4636, explains the previous undetection. We hope that 
systematic deep exposures from ASCA will solve this problem and bring us a complete 
understanding of the mass distribution around elliptical galaxies. 

K.M. acknowledges support by the Postdoctoral Fellowship of the Japan Society for 
Promotion of Science. 
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* The definition of solar iron abundance: Fe/H=3.24xl0 -5 (by number) 

1: Fraction of leakage photons from other annular regions into each ring. 

2: Fraction of stray light, or scattered photons from outside of the field of view into each 

ring. 

Table 1: The fitting results for the ring-cut GIS spectra with the Raymond-Smith model. 



-14- 



Fig. 1 : The 0.7-2.0 keV background-subtracted image of the NGC 4636 field and the south- 
offset field, obtained with the GIS. The contour levels are logarithmically spaced by a factor 
of 1.5, with the lowest contour corresponding to 25% of the brightness of the cosmic X-ray 
background. The image is not corrected for the telescope vignetting. One-dimensional X-ray 
brightness profile, obtained by slicing the image along the dotted region, is also shown. 

Fig. 2: The 0.7-2.0 keV X-ray radial count-rate profile of NGC 4636. Typical /3-model 
profiles convolved with the instrumental response are shown (thin dotted lines). Beyond 
~ 24', the profile shows the data from the south-offset field. We have also shown the profile 
(thick solid line) fitted with a sum of two /5-model components (thick dotted line and thick 
dashed line). The best-fit model parameters are indicated in the figure. 

Fig. 3: Radial mass profiles of the total gravitating matter, the X-ray emitting plasma, and 
the stellar component, in the NGC 4636 system. The solid curves for the total mass and 
the plasma mass assume a constant temperature (0.80 keV) and a constant abundance (1.0 
solar), while the dotted curves take into account the radial gradients in temperature and 
abundance, both for a double /3-model. The dotted region shows acceptable mass curves 
based on the Monte-Carlo simulation. 
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